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ABSTRACT 

The purpose of this Monte Carlo study is to determine the 
accuracies to which the scalar orbital parameters speed, height 
and flight path angle can be determined from 60 seconds of ship 
tracking data. The Apollo Go, No-Go decisionis based upon these 
scalar parameters. Comparisons are included with results ob- 
tained by other investigators using different methods. This study 
shows that speed is very sensitive to the weighting employed, and 
that speed and flight path angle are sensitive to tracking geometry. 
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X~507-67-61 


MONTE CARLO SIMUI.ATION FOR APOLLO 
GO, NO-GO VERIFICATION 

Bernard Kaufman 


I. INTRODUCTION 

When the Apollo spacecraft is inserted into a pai'king orbit about the earth, it becomes im- 
portant to determine the accuracy one might expect in computing the orbit using only a short 
tracking interval from a single shipboard C-band radar. The tracking results obtained by the 
ship are used for verification of spacecraft on boa.rd data received by telemetry for a go, no-go 
decision. Therefore the use of shipboard radar is that of a backup or veriiication system since 
the major decision will be based primarily on the on board position and velocity measurements 
(Reference 1). During the powered flight phase, coverage from ground stations is adequate but 
land based tracking coverage terminates at or prior to cutoff of the S-IV-B stage at about 1440 
nm down-range. The insertion ship thus becomes the primai^y observational base for the insertion 
phase. 

The ”orbit determinations'Lare simulated in this report by utilizing least square fits to 
sampled data generated by Monte Carlo techniques. A measurement error model is first assumed 
and is then i:)ertux'bed by a random process. The deviations are then used in a weighted least 
squares method to determine the "best fit" orbital elements. 

Results that were obtained previously (reference 1) and results not yet published (Table 3) 
were at variance with one another and thus indicated that further work was necessary. This re- 
port was then undertaken as an attempt to resolve these difficulties. 


II. ANALYSIS 

A comment on notation should be made here. Bold face capitals denote matrices while small 
letters denote the elements of the matrix. 

We assume a ra.dar which measures range (r), azimuth (a) and elevation (e) with correspond- 
ing errors S r, 8a and 8e. We may express these errors in terms of the Cartesian position vector 


^ 3 * 1 ) 



i = L2,3 


by the first order terms of the Taylor series 


r + 8r = r + 


3 



+ V, 


1 



where is a residual error term resulting from neglecting higher order terms in the Taj^lor’s 
expansion. 


Then 


Sr - 



+ V, 




Similarly 


y where the summation is assumed 


It must be pointed out here that S x is a time dependent matrix and therefore is not referenced 
to a fixed point in time. But for a least squares fit such a reference point is needed. Now 


^ _ BX 

■ Br 


S r + 


ii 

Ba 


5 a + 


B £ 


S e 


where 


3 X is a 3 X 1 matrix . 


8X 




Sx 

By 

2 


since x,y ,z are not functions of r, a, and e. 
Then we may write 


SX 


(3^1) 




B r 

r 3(x,y, z)' 


Ba 

La(r,a, e) . 

i(3x3) 

_Be_ 


(3x1) 


( 1 ) 


where we consider h\ as being expressed in an inertial coordinate system whose origin is at the 
center of the earth. 

This equation is an inconvenient one to solve in its present form since it is time dependent 
and v/e therefore seek a variational equation in a different form which will allow us to evaluate the 
variation in position at any given time as a function of the measurable s r , a and € . 
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If we define a new coordinate system (z) centered at tlie observer where 


Zj is directed towards the East 
Zj is directed towards the North 
Zj is normal to the local tangent lolane 

then the coordinates of a point in space measured by the radar are given by: 





r s i n a c o s e 


^(3-1) - 

^2 

- 

r cos a cos c 



^3 


r s i n e 






(3x1) 


For the first order terms, we have 





COS e s i n a 

- COS a 

-sin e sin a 


h r 



z: 

cos e COS a 

s in a 

-sin e. cos a 


- r cos € ha 


_*^3_ 


s in € 

0 

cos e 

(3x3) 

The 


or in abbreviated form 

j (3«3) ” ®(3«3) ^ ^(3« 1) 

where J is the orthogonal 3x3 matrix; 


( 2 ) 


( 3 ) 


and 


(3«3) 


1 

0 

0 


0 

- r cos e 
0 


0 

0 


r 


-J(3x3 ) 


SK 


(3x1) 


S r 
ha 


-^(3x1) 


represents noise in the measurements (Reference 2). 

If we let and be the geodetic latitude and longitude respectively of the station then we may 
transform equation (3) to inertial coordinates by means of: 


^(3«l) *^(3«3) (S-S) (^c'> ^(3«l) ~ ^(3’<3) ®(3«1 ) 


3 



where 


\{ 


(3-3) 


(3-3) 



(3-3) 



S( 3 ^i) is the position \^ector of the station in an earth centered coordinate system where Sj is 
tc^wards Greenwich, S 2 is 90"* east ofsj and is along the earth^s axis of rotation; and R^ are 
rotations about the corresponding x and z axis and 9 ^ is the GHA at the time of observation. Us- 
ing this notation equation (3) becomes 


*^3-3) ^^^3-1) "J(^3-3) ^(3-3) ^2t(3-3) ^^(3>1) ^3- 3 ) ^^ ( 3- 3 ) ^ ^(3- 


U 


where SX is inertial. If we include in the above equation the variation of the measurement (Refer- 
ence 3 ) with respect to bias in the measurement (^^^( 3 . 1 )) we may write the equation as 


^(3-3) ^ ^ (3- 1) “ ^(3-3) ^ ^(3- 1 ) J (3-3) ^*(3x3) ^ ^(3- 1 ) “ ^(3-3) ^ ^(3- 1) 


( 4 ) 


where 


^(3-3) ■■ J(3-3) ^(3-3) ^z(3-3) 

Equation (4) permits the computation of for each point in time which is then substituted 

into equation ( 1 ). 

Consider now the position vector X^ 3 xd time T - t, ^( 3 ^ 1 ) related to the state vector 
1 ) ^ functional relationship 

X(3.n "^(3x1) (^0 ) 

(3-U (3-1) ^(3x1) 


and for the erros S\., we have 


X,3.n + S X,.,,, =X.3^,, (X„ +SX(, .Xo +SXq ) 

(3’‘1) (3-1) (3x1) ^(3x1) ^(3-1) ^(3-1) (3-1) 


where we shall consider X.^ to be the state at some fixed ''reference point" in time. Expp.nd- 

0(6-1) 

ing the above as a Taylor series we have 


^ ( 3 - 1 ) ^ ^ (3- 1 ) ^ ^ (3- 1 ) (^0’ ^0) ^ 
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(3- 3) 


"(3-1) 
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SX. 


+ V 


"(3*1) 


(3* 1) 


-O 
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or 


h\ 


(3-1) 


^(x, y. 2 ) 
L^^^'o-yo' ^o-^o-yo' ^o)J 




(6x1 ) 


+ v: 


(3-1) 


(3-6) 


where we do not include the partials of x, y, z since there is no direct measurement of them and 
thus , is the matrix of position only and is the variation in the entire state at the 

reference pWnt. . 

Letting i,. denote the matrix of partial derivatives we have 

(3x6) 


^( 3 * 1 ) ” ^ ^( 3 - 1 ) ^0 




(3x6) 


"(6x1 ) 


or 


(3-1) 


^(3-3) ^ ^(3x 1) 


^0 ^ ^ C 

^^(3x6) ^(6x1) 


( 5 ) 


where is defined as above and ~ L '( 3 xi)** • method of weighted 

least squares with equation (5) where we seek tlie best correction to the state represented by the 
matrix and where L( 3 x 3 >^X^ 3 , n is given by equation (4). The method of least squares 

for the solution of the unknown parameters ^X^ stipulates that: 0 = ^[ 1 . 3 ) ^'( 3 xi) ^ 

minimum 


where 



0 



cos-^ 

0 


0 

p2 rr 2 


(3-3) 


( 6 ) 


is the weighting matrix. 


*This matrix of partial derivatives is colled the state transition matrix when the complete matrix is written as 


B(x, y, z, X, y, z) 

?(^o* yo' ^^ 0 ' yo- ^o>- 


(6x6) 


The computation of this matrix is found in reference 4. 

**lt should be noted that equation (4) may now be written as: 

*"(3«3) 1)' ”3) 1) J (3«3) (3 > 3 ) *®(3 < 1 ) " ^(3 « 3) *P(3« 1) 


where the fixed reference point (^>X,. ) necessary for the least squares procedures. 

(6x1) 


5 


Substituting from (4) 


1 


then 


or 


<1> - [L S X - Lo^ 6 Xp]^ 'V*i [L i> X - LOji « Xq] 

= [SXTIJ _ JXjoJLT] ,v-i [LSX-Lo„ 5XJ 

:: SXT LT W-> L5X - S\l o'f IJ W-' LSX - i XT 17 W'l SX^ + S Xj $t ij y-i 


B<I> 

a(8X„) 


n 0 =. - C.T IT IV-) LS X + 17 W-1 1,0^ 




0 TlT\V-ilSX = «T 17 iv-i 5 Xo 



If we write the above in summation notation where we sum over n observations we have 


2](0TLTffi->L5X)(,.,, 


2[^(0TLTW-i 




(6x1 ) 


or 




''(6x1) 


y (oTLTfi-i Lv 


I — I 
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(6X 6) 




U W-' L5X) 


(7) 


(6x1) 


which are the normalized equations. 

Equation 7 represents the weighted least squares procedure for determining the ’’best fit” 
variations at a fixed reference point. The nonweighted procedure is identical and if one sets 
W " I, the identity matrix, equation (7) reduces to the nonweighted case: 


sx. 


(6x1) 


n 

2](oTO(,) 


y (®j s X) 


J(6x 6) 


Z I 

1 


J(6x 1) 


( 8 ) 


III. TRANSFORMATION TO SPHERICAL ELEMENTS 

Once 5 Xq ^ is obtained either by equation (7) or equation (8) it may be more instructive to 
look at the resulting uncertainties in terms of deviations in the insertion elements of the space- 
craft. Accordingly, we define these elements as: 

h. spacecraft altitude above the earth 
V. speed 

r . flight path angle 
a. insertion azimuth 

I . 

K. right ascension (or longitude) 

8. declination (or latitude) 
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The transformation from injection elements to position and velocity coordinates are as 
follows (References 2 and 5 ); 
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r cos S. cos K. 



1 1 
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=; 

r cos S. s in \. 



1 1 

_z _ 

(3-1) 

_ r s j n S. 


(3-1) 
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y 


(3-1) 


V. (sin 7 . cos S. cos - cosy, sin a. sin - cosy, cos a. sin S. cos\^) 

V. (sin 7 . cos S. sin^. + cos 7 . sin a. cos - cosy, cos a. sin S. sinX. ) 

V. (sin 7 . sin S. + cos 7 . cos a. cos S. ) 


(3-1) 


It can easily be seen that we may then write the transformation for the deviations as 


sx. 


^ (x, y, z, X, y, z) 
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'*'(6-6) P(6-l) 

where is the matrix of partials and is the variational matrix. Then we have 

^( 6 - 1 ) “ '*'(61 6 ) ^ ^ 0 ^ 

and from equation (7) 
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( 11 ) 


Recall from equation 4 that 


I' (3, 3) ^ ^(3, I) 1*(3,3 ) ^ *^(3«l) J (3*3) ^^(3«3) * '^(3» 1 ) 


1 ^( 3 « 3 ) ‘’^( 3 , 1 ) 


If we apply random numbers to the error terms and we then have a means 

of calculating the statistics in a Monte Carlo sense. This is done in the following way; random 
numbers are genei’ated and used as multipliers for the S s 3, and S p 3."^. errors. These-"’ 
products then give modified 28^3,,^ and which are constants over the n observations. How- 

ever for each of the n observations different random numbers are generated for the 5 K , or 
noise errors. Then equations ( 5 ) and ( 7 ) allow us to compute the variations. This procedure 
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is repeated for j times or Monte Carlo samples yielding j matrices for SX^ and p. The standard 
deviations may then be computed. For example 


) 

- (Sh7)2 
j -1 ^ 


1/2 


where Sh- is the mean of the errors in altitude determined from j Monte Carlo runs. Similarly 
for the other injection parameters and for the uncertainties in the state 


IV. CALCULATION OF PERIGEE UNCERTAINTIES 

Deviations in perigee may be included in the Monte Carlo statistics by means of classical 
methods of celestial mechanics. Denote the original ’’unperturbed’^ reference state byX^^. Then 
the ’’best fit” state is 




SX, 


(6x1 


X 

y 


r = (x2 + y2 + z2)l/2 
V = (x2 + y2 ^ 22)1/2 


y 

z 


Then perigee is calculated as follows 



where m is tlie gravitational constant of the earth 


( 12 ) 


e sin E = 



(13) 


which yield the eccentricity e and finally 

rp=a(l-e) (14) 


If perigee is also computed with the original vector then the variation of perigee is thus 
obtained and the standard deviation may be calculated in the usual manner. It must be pointed out 
that for near circular orbits, perigee uncertainties are not gaussian. Thus the standard deviation 
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in this case is not really v/hat it implies. However, once one realizes this, the vStandard deviation 
and the mean for the perigee errors still allow one to intuitively gain insight into what is happen- 


V. EXAMPLES AND RESULTS 

Two different sets of data were selected for numerical examples in a preliminary study of the 
effects of geometry on the uncertainties. Range safety requirements at Cape Kennedy enforce 
certain restrictions on launching and to meet these requirements two bands of 26“^ launch azimuth 
widths were selected. The azimuth spread of 26"^ represents the maximum daily launch window 
for the Apollo mission. 

Two ship locations are necessary to provide coverage of the two azimuth bands. Ship A with 
geodetic coordinates 26° .0 (latitude North) and 47° .5 (longitude West) is used to cover the azi- 
muth band from 72° to 98°. Ship B with coordinates 21° .25 N and 48° .75 W is used to cover the 
azimuth band from 82° to 108°. 

Figure 1 shows the locations of the ships and the two bands of azimuth. The coverage of the 
radar indicated by the circles are for elevation angles above 5 °. 

Table 1 shows the uncertainties for the weighted least squares for height, speed, flight path 
angle and perigee. The last three columns are the average perigee uncertainties and the uncer- 
tainties in total position and velocity. The first 3 columns are launch azimuth, elevation angle of 
the first observation and elevation angle of the last observation. As can easily be seen, the un- 
certainties for Ship A and Ship B are almost identical. The time span for the observations is ex- 
actly one minute along the orbit with the first observation occurring at insertion into the orbit and 
one observation per second thereafter. 

Table 2 is the same as Table 1 except that here the weighting matrix is the identity matrix. 

This has the effect of giving a heavier weight to the angular measurements than is given the more 
accurate range measurement and one would suspect that the resulting uncertainties would be larger. 
A comparison of Tables 1 and 2 shows that the uncertainties in hi, 7 1 and total position are the 
same; however, for speed, the difference is considerable. The deviations associated with perigee 
and total velocity are also lai'ger in the unweighted case. This clearly shows the penalty one pays 
for using nonweighted least squares. 

Table 3 is the results obtained by W. D. Kahn in a study (yet to be published— Ref . 8 ) using a 
linear error analysis program. A comparison with Table 1 for ship B shows that speed, flight 
path angle and the total velocity are in fairly good agreement. However altitude and total position 
do differ by a significant amount. 


P. G. Brumberg in reference 1 shows results that he obtained using Monte Carlo techniques. 

A complete description of the error model used can be found in the reference. A comparison be- 
tween the results obtained by Kahn, Brumbei’g and’the author is shown in Table 4 for the 108 de- 
gree launch azimuth only. The results obtained for Brumberg were read from graphs in reference 
1. These results were obtained without weighting the measurements. 

The above mentioned results show that while the author and Mr. Kahn’s results agree except 
in altitude and position further investigation is still required into this difference. 

As pointed out earlier, the uncertainties in perigee for a near circular orbit are not dis- 
tributed in a gaussian manner. However in the above mentioned examples, the uncertainties in the 
state vector were normally distributed with a zero mean within the expected accuracy of a finite 
sample size. 
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Several other studies he^ve been made in connection with the perigee uncertainties. Reference 
6 utilizes an analytical expression of the go, no-go criterion based on lifetime and results are 
given in that report. A prelimiimry analysis, currently being prepared for publication (Reference 
7), has also been undertaken for the probability distribution of perigee uncertainties and its ap- 
plication to the go, no-go decision. This analysis utilizes two-body equations of motion. 

Figures 2 through 7 show the uncertainties in In , Vi, rp, total position and total velocity 
for Ships A and B plotted versus launch azimuth. Both weighted and unweighted cases were shown. 
The effect of geometry is clearly seen in these graphs. 

The error model used as input in the above examples was as follows: 

Noise (oK) Bias 


Range (meters) 

10.0 

20.0 

Azimuth (milliradians) 

0.4 

0.8 

Elevation (milliradians) 

0.4 

0.8 

STATION ERRORS 

(meters) 


X 

430 


Y 

430 


Z 

0 



WEIGHTING MODEL FOR WEIGHTED CASES 
a range = 10.0 meters 
a azimuth = 0.4 milliradians 
cr elevation = 0.4 milliradians 
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75 
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.90 
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.026 
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1.23 

6.39 
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Bios (i‘p) Noise (;K) Station errors {meters} 
r (meters) 20.0 10.0 x 430.0 

a (millirodions) 0.8 0.4 y 430.0 

e (milliradla.ns) 0.8 0.4 2 0.0 
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-4.7 

1.24 

6.79 

Ship B 

82 

7.6 

10.1 

.94 

2.70 

.024 

5.1^ 

-4.3 

1,28 

6.76 

85 

8.G 

13.9 

.90 

2.23 

.02 9 

4.1 

-4.0 

1.23 

7.15 

90 

11.8 

25.7 

.80 

1.48 

.038 

3.2 

-4.0 

1.09 

7.80 

95 

13.2 

40.2 

.77 

.731 

.046 

3.4 

-4.5 

1.04 

8.35 

97 

13.0 

36.2 

.78 

.890 

.044 

3.3 

-4.3 

1.04 

8.10 

100 

11.7 

25.4 

.81 

1.50 

.037 

3.3 

-4.1 

1.07 

7.50 

105 

8.5 

13.7 

.92 

2.29 

.029 

3.8 

-4.4 

1.20 

6.93 

108 

7.5 

9.8 

.96 

2.78 

.023 

4.5 

-4.8 

1.25 

6.65 

Bias (5p) Noise (5K) Stolion errors (nerers) 
r (meters) 20.0 10.0 x 430.0 

a (n’.i 1 1 irac'ian s) 0.8 0.4 y 430.0 

£ (ml 1 ! irodions) 0.8 0.4 z 0.0 
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Tablo 3 

Linear Error rViin lysis (ieeni) 


OiK.- sinf.io lince; nii.’i ! ! :: s f:-; Apcjllo Go, Ko-G^> v'-rifU ntioi’ 
AltlTUCC loSA.^ !xni 
Gcco/,;i it i i/ .00 1 io .00/ 

Loclvioy fioni ir;'//!;or; {/ orio nOnuiO (or fro/i 

CJ r/ininMin; cl'; /o 1 i on of 7.5' 1o c>no 'Tiinuto !cMor) 

Ship A 26 LO (Gcccl. I.ot, Gorih); 47 .5 (Lone. V/c-st) 
Ship B 2LL25 (Gooa. Lee Norih); 48'' .75 (Long. V/esl) 


LaiMch 

A/.irmilii 



1 

' 1 

Total 

lAisition 

(l^m) 

To La] 
Vole city 
(m/.s) 



SI lip D 



82 

. () 3 

.00 

.024 

1.11 

6.85 

85 

.[>0 

.50 

.025 

.95 

6.91 

50 

.30 

.41 

.030 

.76 

7.33 

95 

.18 

.41 

.040 

.68 

8.25 

100 

.30 

.42 

.030 

.76 

7.33 

108 

.63 

.G7 

.024 

1.11 

6.87 



Bios (bp) 

Noise (oK) 

Stc-ion cnors (meters 

r (meters) 

20. 0 

10.0 

X 430.0 

a (mi I lira d inns) 

0.8 

0.4 

y 430.0 

e (mi 1 1 irodi ons) 

0.8 

0.4 

2 0.0 


Table 4 

108 Launch Asiinuth 

One sigma uncertainties for Apollo Go, No-Go verification 
Altitude 185.2 Km 
Eccentricity .001 to .002 

Tracking from insertion to insertion plus or-.e minute (or from 
c rMiriimurn Cicvction angle of 7.5 to one rriinute Ictor) 
Ship A 26"" .0 (Geod. Lot. North); 47 '‘.5 {Long. V/est) 
Siiip B 21 .25 (Good. Lot. North); 48 .75 (Lor^g. V/eot) 



" ~1 

(ni/scc) 


(len) 

Total 

Position 

(Icni) 

Total 

V'olocUy 

(in/sec) 

FCahn 

.67 

.024 

.63 

1.11 

6.87 

}3ruml;org 

2.1 

.03 

.63 

1.5 

8.3 

Ktiuiman 

(weighted) 

.59 

.023 

.96 

1.25 

6.04 

Kaufman 

(unvveigiiicti) 

2.78 

.025 

.96 

1.25 

6.G5 

L._„ 


LATITUDE (degrees) 












ONE SIGMA HEIGHT UNCERTAINTIES (Km) 







ONE SIGMA SPEED UNCERTAINTIES ( 






. 0:>0 



o 

lU 

o 


.01 



One Sfg.-ia uncer fa m fie s for 

AIttfucJe 1S5.2 Krr= 

Eccenfricify . 00 ) fo .002 


Apollo Go, P<c-Go vcrfficction 


TfO:<,ng from inserl.cn fo insertion p-fu 
o r.’tr.irnun elevofion angle of 7.5’ 

SMpA 26 '.C (Gcc-J. Lct.^NoTTV 47 
SH.p B 2)'. 25 (Gccd. Laf. Norifo' 43' 


one fTunufe for from 
to one minute Icfcr ) 
.5 fLcng, V.'esfi 
■75 iLcng V.'cs*i 


^ BiASiSil 

20.0 

Qfrnrcdl 0.8 

<fmrc4) 0 8 

STAT.L^ (m. 

> <40.0 

y <30.0 

2 0.0 


^OISE (3 5) 
fO.O 

0 4 

0.4 




J... 

75 


80 



launch azimuth (deg) 

Figure 4 




ONE SIGMA PERIGEE UNCERTAINTIES (Km) 



1.5 


1.0 


0.5 


One sigrno unccrfointies for Apollo Go, hiO-Go verificcjic 
Altitudo 185.2 Km 
Eccentricity .001 to .002 

TracUng from insertion to insertion pi'js one minute (or fi 
G riinimum clcvojicn angle of 7 5^ to one rr.ii'.jte la 

Ship A 26 .0 (Good. Let, Ncrrh), 47 .5 (Long West 
Ship B 21\25 (Good. I.ot, Ko"th), 43=^ .75 (Long. Y/est 


r (.-n) 

Oi.mrccf) 

<<mrcdi 


BIAS KC'tsr (i^) 


20.0 

0.0 

o.e 


10 . c 
o.< 
o.'i 


STATiWt Efx-HORS tm) 
X 4i0.0 

y 430.0 

2 0.0 
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Figure 5 
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ONE SIGMA TOTAL POSITION UNCERTAINTIES (Km) 



Figure 6 
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TOTAL VELOCITY UNCERTAINTIES 







